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ABSTRACT 

We present a detailed comparison between nu merical cosmological hy drodynamic 
zoom simulations and the semi-analytic model of lSomerville et al.l ()2008a| ). run within 
merger trees extracted from the simulations. The high-resolution simulations repre- 
sent 48 individual halos with virial masses in the range 2.4 x 10^"'^ Mq < Mnaio < 
3.3 X IO^'^A/q. They include radiative H & He cooling, photo-ionization, star forma- 
tion and thermal SN feedback. We compare with different SAM versions including 
only this complement of physical processes, and also ones including supernova driven 
winds, metal cooling, and feedback from Active Galactic Nuclei (AGN). Our analysis 
is focused on the cosmic evolution of the baryon content in galaxies and its division 
into various components (stars, cold gas, and hot gas), as well as how those galaxies 
acquired their gas and their stellar mass. Both the SAMs and simulations are com- 
pared with observational relations between halo mass and stellar mass, and between 
stellar mass and star formation rate, at low and high redshift. We find some points 
of agreement and some important disagreements. SAMs that include the same physi- 
cal processes as the simulations reproduce the total baryon fraction in halos and the 
fraction of cold gas plus stars in the central galaxy to better than 20%. However, 
the simulations turn out to have much higher star formation efficiencies (by about a 
factor of ten) than the SAMs, despite nominally being both normalized to the same 
empirical Kennicutt relation at z = 0. Therefore the cold gas is consumed much more 
rapidly in the simulations and stars form much earlier. Also, simulations show a tran- 
sition between stellar mass growth that is dominated by in situ formation of stars to 
growth that is predominantly through accretion of stars formed in external galaxies. 
In SAMs, stellar growth is always dominated by in situ star formation, because they 
significantly underpredict the fraction of mass growth from accreted stars relative to 
the simulations. In addition, SAMs overpredict the overall gas accretion rates relative 
to the simulations, and overestimate the fraction of "hot" relative to "cold" accretion. 
We discuss the reasons for these discrepancies, and identify several physical processes 
that are missing in our SAM and in other semi-analytic models but which should be 
included. We also highlight physical processes that are neglected in the simulations 
studied here, but which appear to be crucial in order to understand the properties of 
real galaxies. 

Key words: galaxies: formation; galaxies: evolution; methods: N-body-simulations; 
methods: numerical 
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The dark energy dominated dark matter paradigm (ACDM) 
provides a successful theoretical model for understanding 
and simulating galaxy formation. Within this framework the 
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large-scale structure of dark matter develops from initially 
small-scale density fluctuations, i n a bottom-up scenari o 
only driven by gravitational forces ijBlumenthal et al]ll984h . 
The baryonic component follows the dark m atter, is shock- 
heated, cools and condenses into galaxies l|White fc ReesI 
Il978al :[ Rees fc Ostriker|[l977l ). 

The assembly of dark matter halos, which domi- 
nate the total matter content in the Universe, in large 
cosmological volumes can be followed with merger trees 
based on analytic approaches, e.g. using Monte-Carlo 
methods based on the extended Prcss-Schcchter formal- 
ism (EPS, Press & Schochtcr 1974; Bower 1991 ; Bond ct al. 
199ll; ISomerville fc Kolatt] Il999l: [Neistein fc Dekell l2008l : 



Zhang et all |2008| ; lAngulo fc Whitd bOloD . Alternativelv. 
the full dynamical evolution of dark matter can be ac- 
curately followed with coUisionless particles in direct nu- 
merical simulations which are, by now, well resolved at 
the relevant scales (|Frenk et al ] 1 19881 ; iNavarro et al. 1 1 19971: 
Moore et af] Il999l; iKlvpin et al.l Il999l; [ Bode fc Ostrikerl 



20031; ISprineel et al.l |2005| . |2008| ; iDiemand et all l2008l : 
Klvpin et al.ll2010t l. Here the identification of dark matter 
halos and the construction of merger trees is considerably 
more demandin g and various differ e nt approaches have been 



discussed (e.g. I Davis et al!" 'l985'; 'Kauffmann et al.' ' 19931; 



Ghigna et al.1 12000; Springcl ct al. 2001; WcUcr ct al. 200l 
_Genel et al.''2008; Fak houri fc Mai 120081 ; iPlanelles fc Qiiilij 
201Q; iSkorv et al. 2010l a nd references therein). In particu- 
lar, Springel et al.l l|2005l ) constructed and analysed merger 
trees for the Millennium s imulation, using t he "Friends of 
Friends" (FOF) t echnique (jPavis et al.ll 19851 ') to identify ha- 
los, and SUBFIND (|Springel et al. !'2001'l to identify sub-halos 
(bound objects within larger virializod dark matter halos). 

Simulations of the formation and evolution of the galax- 
ies which are believed to inhabit these dark matter ha- 
los are more demanding, theoretically as well as numer- 
ically. Additional gas-dynamical and radiative processes, 
such as the formation of stars and black holes as well 
as the respective feedback, have to be taken into ac- 
count. To follow the evolution of galaxies two main ap- 
proaches have been developed over the past decades: Semi- 
analytic mod els fSAMs) and direct cosmological simula- 
tions. SAM s IWhite fc Frenk" 'l99l'; 'Kauffmami et al.' ' 19931; 
Cole et al.1 [l9 94: Kauffmann 1996; Somcrvillc fc- Primack 
19991: iKa uffmann ct al. 1999; Kauffmann fc Hachnclt 2000; 
Cole et al. 2 000: Sprin gel et al. 2001; Hatton et al ,2003 ; 
Kang et al.|[20 05: Bau gh et al.ll2005l:lKhochfar fc Silkll2006al; 
Croton et al. 2006; Bower et al.' '2006'; 'De Lucia fc Blaizot 
2007; Somcrvillc ct al. 2008a; Font ct al. 2008; Guo ct al. 
20091 : IWeinmann et al.l 120091 ') use pre-calculated dark mat- 
ter merger trees either from EPS or direct cosmological 
simulations and follow the formation of galaxies with sim- 
plified, physically and observationally motivated, analytic 
recipes. The computational costs of this approach are typ- 
ically low, and the influence of different physical mecha- 
nisms can be investigated separately in a straightforward 
way. Modern SAMs are quite successful at reproducing ob- 
served statistical properties of galaxies in large cosmologi- 
cal volum es over a large range of galaxy masses and red- 
shifts fe.g. ISomerviUe et al.|[2008al : ICuo et al.ll2010l '). Disad- 
vantages are that the dynamics of the baryonic component 
(gas and stars) and the interaction between baryonic matter 
and dark matter are not followed directly and that in many 



cases the assumed models are simplified and use a large num- 
ber of free parameters to fit diffe r ent observations simul- 
taneously (ISomerville et al.1 l2008al ; iBenson fc Boweij l201ll : 
iBower et al.ll2010h . 

Direct cosmological galaxy formation simulations can 
follow the evolution of dark matter and gas explicitly. Even 
though they treat the underlying dynamics more correctly 
than SAMs, the spatial and mass resolution, at present, 
is not high enough to accurately simulate intermediate 
and low mass galaxies in large cosmological volumes. In 
addition, small scale processes like e.g. the formation of 
stars and black holes with the associated feedback has to 
be computed in a si mplified manner with su b-grid/sub- 
resolution models ((Cen fc Ostriker 19931; Dave et "al] 



20011; ISpringel fc Hernguistl l2003l; iMaller fc Bullockl |2004|; 



Nagamino ot a l.l |2005|; iKeres et al.l 2005al : Navarro et al.l 
2009; Schavc et al.l 12010a '). which again require the in- 



troduction of parameters. Ab initio cosmological zoom 
simulations with proper cosmological boundary conditions 
enable direct simulations of the baryonic physics of certain 
regions of interest at higher resolution, either limited to 



popular Iv. individual halos (INavarro & Steinmetz 1997; 


Governato et al. 20071: Naab et al.1 |2007|; 


Brooks et al.l 


2009al: lOser et al.l 2010l: IWadepuhl fc Springell l201ll; 


Puchwcin et al.l 2010|; |Tcvssicr et al.| 2010l: 


Sawala et al.l 


2010; Piontek fc Steinmetz 2011;,Agertz et al.l 


20111). These 



simulations can attain very high resolution, and provide a 
way to resolve galaxies of very different masses with the 
appropriate resolution in each case. However they are very 
time consuming and therefore not currently feasible for 
representative studies of large populations of galaxies. In 
addition, the sub-resolution models are uncertain and it 
is still unclear how sensitive various results may be to the 
details of these sub-grid models or the parameter values. 

Both approaches make definite predictions for the evo- 
lution of galaxy properties at various masses over cos- 
mic time. Because of their greater computational efficiency, 
SAMs generally include more models for physical processes 
than current numerical simulations, and because of their 
greater flexibility, it has been possible to tune them to ob- 
tain quite good agreement with a broad range of galaxy 
properties in the local Universe. SAMs have also been 
shown to reproduce the statistical properties (e.g. lumi- 
nosity and stellar mass functions, star formation rates) of 
high redshift galaxies [z ^ 6) quite well, at least for mas- 
sive galaxies (m star r 
iFontanot et al.1 2009) 



, 10"" M©; e.g. ISomerville et al.1 [20T1I ; 

Therefore we might expect SAMs to 



do a better job of reproducing the observed universe than 
the simulations, but we might worry that they could do so 
for the wrong reasons. Because there is a great deal of un- 
certainty in many of the important processes, and most of 
the physical recipes contain free parameters, if one physical 
process (e.g. gas cooling and accretion) is modelled inaccu- 
rately in the SAM, it is currently possible to compensate by 
tuning a competing process (such as feedback). By running 
SAMs within merger trees extracted from numerical hydro- 
dynamic simulations, in order to constrain the evolution of 
the dark matter component to be the same in both cases, 
we can isolate various physical processes and attempt to im- 
prove the accuracy of the semi-analytic recipes. At the same 
time, by comparing the detailed predictions of the forma- 
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tion histories of galaxies in the simulations with the SAM 
predictions, we may gain insights into the origin of existing 
discrepancies between the simulations and the real universe. 

Various comparison studies between simulations and 
SAMs for large galaxy populations as well as individual halos 
have been discussed in the literature (we summarize these re- 
sults in section [2]) following different philosophies. For some 
studies only i ndividual physical processes, like coolin g, were 
investigated (ILu et al.ll201ol : iBenson fc BoweJ 120111 ). while 
others focussed on the ev olution of individ ual objects, as a 
high-m ass galaxy cluster llSaro et al.ll2ofoh or a single disk 
galaxy (|Stringer et"al]|2010l 'l. 

Our approach is new in many respects. We compare 
the evolution of individual halos but use the, up to now, 
largest num ber of high - resolu tion zoom simulations (48), 
presented in lOser et al.l (|2010l ). The simulations cover dark 
matter halos in the mass range of 2.4 x IO^^A^q < A^Haio < 
3.3 X IO^^Mq. Although a limited complement of physical 
processes have been taken into account in the simulations, 
the more massive of these halos have been shown to repre- 
sent fairly well t he evolution of observed massive galaxies 
l|Oser et al.ll2O10D . These simulations are c o mpare d to re- 
sults from the full SAM of ISomerville et all (|2008al '). which 
has been shown to represent present day galaxy properties 
reasonably well over a wide range of masses, as well as dif- 
ferent stripped down versions. The SAMs are run within 
merger trees extracted directly from the numerical simula- 
tions. In addition, we compare both model predictions to 
observations at different redshifts and point out where the 
respective models succeed or fail either to match each other 
and/or the observations. 

The paper is organized as follows: In section[2]we discuss 
results from previous comparisons between SAMs and sim- 
ulations. The hydrodynamical simulations and the merger- 
tree construction method used for this study are discussed 
in section [3] and we br iefly review the ingredients of the 
ISomerville et all (1200831 1 SAM in sectional The redshift evo- 
lution of the baryonic components in simulations and SAMs 
is compared in section [5] followed by a comparison to obser- 
vations in section [S] In section [7] we summarize and discuss 
our main results. A resolution study for individual halos can 
be found in the Appendix. 



2 PREVIOUS COMPARISON STUDIES 

Previous quantitative comparisons between simulations and 
SAMs have either focused on whole populations of galaxies 
llBenson et al ''200l': "V oshida et al."2002': 'HcUv ct al."2003; 
ICattaneo eTal...2007. : .Benson fc Bower.. 201 1: .Lu et al.. .201 
or individual objects like a galaxy cluster ijSaro et aLll201 
and a single disk galaxy ("Strin ger et al.|[2010l V Hellv et al.l 
l|2003h compared the efficiency of gas cooling for differ- 
ent dark matter halos as a function of redshift between a 
(50Mpc//i)^ SPH (Smoothed-Particl e-Hydrodynamics) sim- 
ulation (Hydra, iPearce et al.l I2OOII ) excluding star forma- 
tion, heating and feedback and a stripped down version - 
without star form ation or feedback - of the Galform SAM 
l|Cole et al.ll2000l ). For z — they find good agreement of the 
cold gas mass between the SPH simulation and the SAM. 
At high redshifts, however, more gas tends to cool in low- 
mass halos in the simulation due to the limited numerical 



resolution. Still, they conclude that simulations and SAMs 
give consistent results for the evo lution of cooli r ig gal actic 
gas and confirm earl ier findings of I Benson et al.l (120011 ) and 

lYoshidaet all hOOt . 

ICattaneo et al.l (|2007l ) did consider star formation 
and supernova feedback for a similar comparison in a 
34.19Mpc'^ volume. For t he SPH-simulation they used 
TreeSPH (|Dave et al] 1 19971 and the input for their SAM- 
code GallCS (|Hatton et al.ll2003l ') was the merger trees con- 
structed from the dark matter component of the SPH- 
simulation. The SAM did not include a photo-ionising back- 
ground but followed the cooling by metals, while the simula- 
tions did include a photo-ionizing background and assumed 
primordial composition of He and H. The star formation 
prescription in the SAM was quite standard (star forma- 
tion occurs above a certain gas surface density, according 
to a Kennicutt-Schmidt-like relation), and a simple recipe 
for SN-driven winds was also included. In order to repli- 
cate AGN feedback, star formation is quenched when the 
bulge component of a galaxy reaches a critical mass. For 
the comparison they used two different SAM versions: one 
with no feedback and the 'full' model. In general, they found 
good agreement between simulations and the no-feedback 
model for the baryonic mass functions at different redshifts 
and in different environments. Moreover, simulations and 
the no-feedback SAM made similar predictions for the 'hot' 
and 'cold' mode gas accretion histories of galaxies (e.g. 
iKeres et al.]|2005bh . However, at low redshifts, much less gas 
was left over in the simulation than in the no-feedback SAM 
with both approaches over-predicting the observed baryonic 
mass function, in particular at the high mass end. The full 
SAM, on the other hand, matched the observations due the 
inclusion of supernova-driven outfiows and AGN feedback, 
which suppresses gas cooling in large halos. They concluded 
that the simulations and the no-feedback model failed as a 
consequence of missing physics rather than computational 
inaccuracies. 

ISaro et al] (|2010t ) compared the galaxy populations 
within a massive cluster (A/ciustcr ~ 1.14 x 10^^ Mq) us- 
ing a high- resolution cosmolo gical re-simulation run with 
GADGET2 JPolag et all 120091') and the SAM model of 
iDe Lucia fc Blaizoj ( 2007 '). They focused on differences be- 
tween the central and the satellite galaxies considering only 
gas cooling and star formation and neglecting any form of 
feedback. In general, they find similar statistical properties 
for the galaxy populations, e.g. the stellar mass function 
with a few remarkable object by object differences. The cen- 
tral galaxy in the simulation starts with a more intense and 
shorter initial burst of star formation at high redshift and 
forms fewer stars at low redshift than in the SAM. While 
in the SAM all stars in the central galaxy are formed in 
its progenitors, in the simulations the final stellar mass is 
larger than the sum of all progenitors. Satellite galaxies can 
lose up to 90 per cent of their stellar mass due to tidal 
stripping - a proces s , whic h is, however, not included in the 
iDe Lucia fc BlaizotI (|2007l ) semi-analytic model, nor in most 
models discussed in the recent literature. 

Moreover, IStringer et al.l (|2010t ) presented a compari- 
son for the evol ution of a single disk galaxy using the SPH- 
code Gasoline (|Wadslev et al.l 20041 ') and the semi-analytic 
model Galform ( Bower et al.ll2006l ) based on the dark mat- 
ter merger history of the simulation. They find that the two 
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techniques show a potential consistency for the evolution 
of the stellar and gas components by assuming the same 
physics and the same initial conditions. They try to mimic 
in the SAM the 'blast wave' SN feedback implemented in 
the simulation, i.e. after a supernova explosion no cooling is 
allowed in a certain v olume. However, u sing the Galform 
model as described in lBower et al.l l|2006f ) (including chemi- 
cal enrichment, supernova and AGN feedback), the resulting 
system is not recognisably the same as predicted by simu- 
lations. At all redshifts, the stellar mass is much larger and 
the hot gas fraction is much lower in the simulation than in 

the SAM. 

Finallv. lLu"et al.l l|2010l ) and lBenson fc Boweij l|201lf ) fo- 
cus on the algorithms for ga s cooling in SAMs in great de- 
tail, n enson fc Bowed iioij) compare cold (rapid) and hot 
(slow) accretion rates in the G alform SAM and in simu- 
lations from Keres et al.l l|2009t ) (50 Mpc//i, 2 x 288^ parti- 
cles). They used their 'full' model including feedback and 
metal cooling, although these processes are not included 
in the simulations. Moreover, they modified their SAM by 
adopting an updated calibration for the transition between 
the rapi d and slow cooling regime following the method- 
ology of iBirnboim fc Dekell (|2003ah . They find reasonably 
good agreement for the hot and cold mode accretion frac- 
tion in the SAM and the simulations and thus, they conclude 
that the cold-mode physics is already adequat ely accounted 
for in SAMs. In the stu dv oflLu et al.l (|2010| ') five different 
SAMs ('Munich' model: ICroton et al.ll2006l. ' K ang' m odel: 
iKang et al.ll2005l. 'Galform' in odel: ICole et al.]|2000l . 'Gal- 



ICS' model: 



Hatton et al. [j2003l and the 'Somerville' model: 



ISomerville fc Primack 1999') are compared to the simula- 



tions of iKeres ct a l. (2009), without considering any feed- 
back or metal enrichment in either method. They find a 
significant difference between hot and cold accretion rates: 
compared to the simulations, the cold mode accretion rates 
are too low and the hot mode accretion are too high in 
SAMs. They construct a modified cooling recipe for the SAM 
to enable simultaneous hot and cold accretion, resulting in 
much better agreement between the SAMs and the simula- 
tions. 

Throughout the course of this paper, we will refer back 
to these studies and comment upon the similarities and dif- 
ferences with our results. 



3 THE SIMULATION AND MERGER TREE 
CONSTRUCTION 

3.1 Simulation setup 

The cosmological zoom simula tions presen t ed in this pa- 
per are described in detail in lOser et ^l] l|2010l ') and we 
briefly review the simulation setup here. The daxk mat- 
ter halos for further reflnement were selecte d from a dark 
matter only N-body simulation (Gadget-2, ISpringel et al.l 
I2005D with a comoving peri odic box length of L = 100 Mpc 
and 512^ particles (see also lMoster et allboiol ). We assume 
a ACD M cosmology b ased on the WMAP3 measurements 
(see e.g. ISpergel et all 2003) with as = 0.77, fi™ = 0.26, 
Qa = 0.74, and h = Ho/(100 kms~^) = 0.72. The simula- 
tion was started at 2 = 43 and run to 2: = with a flxed 
comoving softening length of 2.52 h~^kpc and a dark mat- 
ter particle mass of Mdm = 2 x lO^M©/^. Starting at an 



expansion factor of a = 0.06 we constructed halo catalogues 
for 94 snapshots until z — separated by Aa — 0.01 in 
time. From this simulation, we picked 48 halos identified 
with the halo finder algorithm FOF at 2: = 0. To construct 
the high-resolution initial conditions for the re-simulations, 
we trace back in time all particles that are closer than 2-r2oo 
to the center of the halo in any snapshot and replace them 
with dark matter as well as gas particles at higher reso- 
lution (r^i, = 0.044, fioAf = 0.216). In the high resolution 
region the dark matter particles have a mass resolution of 
muM = 2.1 • IO^Mq/i"^, which is 8 times higher than in the 
original simulation, and the gas particle masses are mcas = 
T7^sta^ = 4.2 ■ IO^Mq/i"^. Individual cases were run at 64 
times higher mass resolution and 4 times higher spatial res- 
olution. The re-simulated halos cover a mass range of two or- 
ders of magnitude (2.4 x lO^M© < Muaio < 3.3 x lO^M©). 

For modeling the gas component we u se the entropy 
conserving formulation of SPH (Gadget-2, ISpringel et al.l 
l2005h . We include star formation and cooli ng for a pri- 
mord ial composition of hydrogen and helium IjTheuns et al.l 
Il998l ). The cooling rates are computed under the assumption 
that the gas is optically thin and in ionization equilibrium. 
Furthermore, our simulations include a spatially uniform 
re dshift dependent UV back ground radiation field according 
to lHaardt fc Madaul l| 19961 ). where reionization takes place 
at 2 « 6 and the radiation field peaks at 2 « 2 — 3. 

To mo del star formation and SN f eedback we use the 
approach of lSpringel fc Hernguistl (|2003l ). In this model, the 
ISM is treated as a two-phase medium where clouds of cold 
gas form from cooling of hot gas and are embedded in the 
hot gas phase assuming pressure equilibrium. The hot gas is 
heated by supernovae and can evaporate the cold clouds. 
Stars form from the cold gas whenever the local density 
exceeds a threshold density (p > pth = 0.205cm~^). The 
star formation rate is calculated by 



(1) 



Here, /3 is the mass fraction of massive stars, which are as- 
sumed to explode as supernovae type II, pc is the density 
of cold gas and t* = t1{p/ pth)~^^'^ is the star formation 
time scale. The supernova explosions heat the surrounding 
gas w ith an energy input of lO^'^ ergs. ISpringel fc Hernguistl 
(l2003l ) used an idealized, isolated disk galaxy simulation to 
set the free parameters pth and t^, by adjusting them to ob- 
tain a match to the observed Schmidt-Kennicutt relation. 
We adopt the same values of these parameters here. 



3.2 Merger trees 

We extract the merger trees for the dark matter component 
directly f rom the c osmological re-simulations as described in 
iHirschma nn et al. 112010). For every snapshot at a given red- 
shift, we first identify individual dark matter haloes using 
a FOF (Friends-of-Fri ends) algorithm w ith a linking length 
of 6 = 0.2 (fa 28kpc, iDavis et al.lll985h . In a second step 
we extract the subhalos of every FOF group using the Sub- 
find algorithm (Springcl ct al. 2001). This halofinder identi- 
fies over-dense regions and removes gravitationally unbound 
particles. In this way we split the FOF group into a main or 
host halo and its satellite halos. In most cases, 90% of the 
total mass is located in the main halo. 
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Figure 1. Visualisation of merger trees for four re-simulated halos with different masses: upper leftiMvir = 8 X IQI^Mq, upper right: 
Mvir = 1 X IO^^Mq, lower left: Mvir = 5 X IO^^Mq, lower right: Mvir = 1 X IO^^Mq. Black circles show the dark matter halo at every 
time-step of the simulations. The symbol size is proportional to the square root of the halo mass normalized to the halo mass at z=0. 
The yellow stars indicate the stellar mass, the blue and red filled circles the cold and hot gas mass within the virial radius of the dark 
matter halo. The symbol sizes for the baryons scale with the square root of the masses normalized to the maximum total baryonic mass 
at z=0. 



The sizes and virial masses of the main halos (i.e. the 
most massive Subfind halos) are determined by a spherical 
overdensity criterion. The minimum halo mass is set to 20 
particles (5x10* Mq /h). In the following, we will use isolated 
merger trees which are constructed only for the the main 
halos, i.e. the central objects of one FOF group identified by 
Subfind. The mass of a central object is defined by the dark 



matter mass within the virial radius using the overdensity 
approximation in the sp herical collapse model according to 
iBrvan fc NormanI (|l99a) . The algorithm to connect the dark 
matter halos betwee n the snapshots at differ ent redshifts is 
described in detail in lMaulbetsch et al.l |20o3). The branches 
of the trees for 2 = halos are constructed by connecting the 
halos to their most massive progenitors (MMP) at previous 
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snapshots. Thereby, halo j with Uj particles at redshift Zj 
with the maximum probability p(i, j) is chosen to be a MMP 
of halo i containing rii particles at redshift Zi (where j < i). 
The probablity p{i,j) is defined as 



nov{i,j) 



with 



(2) 



riov = ni{zi) nnj{zj) and 
nmax{i,j) — max{ni{zi),nj{zj)) 

Here, the number of particles found in both halos 

and Umax is the particle number of the larger halo. We re- 
move 'fake' haloes which exist only within one timestep and 
have no connection to any branch (halo masses are gener- 
ally near to the resolution limit). The low redshift ends of 
the branches are then checked for mergers. A halo j at tj 
is assumed to merge into halo i at ti, if at least 50% of the 
particles of halo j are found in halo i. In case of a merger 
the branches are connected. 

Note that the tree-algorithm is only applied to the dark 
matter particles - star or gas particles are not separately 
traced back in time. They are assumed to follow the evolu- 
tion of the dark matter. Therefore, we assign to each dark 
matter halo in a tree a hot /cold phase gas mass by counting 
hot/cold gas particles within the virial radius of the cen- 
tral halo. The stellar and cold gas particles within 1/10 of 
the virial radius are defined as the stellar and gas mass of 
the central galaxy. We distinguish between hot and cold gas 
particles by using the following definition: 



logT < 0.31ogp + 3.2 -^cold 
logT > 0.31ogp-f3.2 ^hot 



(3) 
(4) 



Note that the above discrimination between hot and cold gas 
was established by looking directly at the phase diagrams of 
the re-simulations, where we have divided between the gas 
in the disk heated by SN feedback and the shock heated gas. 
With the above definition for cold gas we mainly capture the 
dense, star-forming gas. 

In Fig. [1] we show a visualization of four merger trees 
of re-simulated halos with virial masses of 8 X 10^2 Mq, 
1 X IO^^Mq, 5 X 10" Mq and 1 x lO^Mg. The size of the 
black circles approximates the dark matter halo mass, the 
yellow stars the stellar mass within the virial radius and the 
blue and red filled circles the cold and hot gas component, 
respectively. The symbol sizes scale with the square root of 
mass normalized to the final dark matter halo mass (dark 
matter component) and to the final baryonic mass (star, hot 
and cold gas mass). We clearly see that galaxies at high red- 
shift contain more cold gas, which either turns into stars 
or is heated towards lower redshifts. In general, for more 
massive halos the fraction of cold gas and stars at 2: = is 
lower. 

To study the influence of numerical resolution on the 
evolution of the dark matter and the baryonic components 
we have simulated a few halos with 4 x higher spatial resolu- 
tion (= 64 X higher mass resolution) than the original dark 
matter simulation. A comparison of the results can be found 
in the Appendix. The overall mass assembly of the main ha- 
los and the number of major mergers do not show any sig- 
nificant variation, although the number of identified minor 
mergers increases due to the higher resolution. Overall, we 



conclude that our results are well-converged and would not 
change significantly if we improved the resolution. 



4 THE SEMI-ANALYTIC MODEL 

The merger-trees constructed as described above are 
used as input f o r the semi-analytic model described in 
ISomerville et all (j2008al . hereafter SOS). The SAM makes 
use of merger trees for "isolated" halos only, and treats the 
evolution of sub-structure within virialized halos using semi- 
analytic approximations. The 'full' version includes photo- 
ioniziation, gas cooling, star formation, SN feedback, metal 
enrichment, and black hole growth in a radio and quasar 
mode with corresponding feedback. However, to provide a 
more meaningful comparison to our simulations, we do not 
only consider the 'full' version, but also 'stripped down' 
models by separately switching off AGN feedback, metal 
cooling. Supernova-driven winds, and 'thermal' Supernova 
feedback. We consider the following different versions: 

• NF: no Feedback, primordial metallicity 

• SN: thermal SN-feedback, primordial metallicity 

• SNWM: thermal SN-feedback, SN-driven Winds, 
metal cooling 

• FULL: 'full' version, including thermal SN-feedback, 
SN winds, metal cooling, and AGN feedback 

In the following we brieffy summarize how the different 
physical mechanisms are implemented and how they differ 
from the ones in the simulations. F or full details we refer the 
reader to lSomerville etal] (|2008al ). In Table [1] we provide a 
summary of the galaxy formation parameters used here. In 
Table [2I we give an explicit overview of the physical recipes 
assumed in the different SAM versions starting from the 
NF model. These are compared to the physics which are 
implemented in simulations. 

(i) Radiative cooling: The rate of gas condensation via 
atomic cooling is comp uted based on the model proposed by 



atomic cooling is comp uted based on the model proposed 
IWhite fc FrenkI l|l99ll ). The cooling time is computed as 



tc 



3/2^mpkT 
p,{r)A{T,Zn) 



(5) 



Here, T is the virial temperature, /imp is the mean molec- 
ular mass, Pg(r) is the radial density profile of the gas and 
A(r, Zh) is the cooling function, which is temperature and 
metallicity dependent. The cooling time is the time required 
for the gas to radiate away all its energy starting at the virial 
temperature. The gas density profile pg{r) is assumed to fol- 
low an isothermal sphere: Pg{r) = mhot/(47rri,ir7'^). Putting 
this expression in Eq. [5] one can solve for a cooling radius 
f'cooi- Within the cooling radius all gas can cool within the 
cooling time tcooi- The cooling rate for the mass within rcooi 
is 



dt 



1 rcooi 1 
^jTOhot :: . 

^ ^vir ^-cool 



(6) 



Following ISpringel et all (|200ll ) and ICroton et aP (|2006D 

it is assumed that the cooling time is equal to the halo dy- 
namical time tcooi = idyn = fvir/Vvir. Two different modes 
of accretion are distinguished: the rapid ("cold mode") and 
the slow ("hot mode") cooling regime. In the rapid cooling 
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Table 1. Summary of the galaxy formation parameters in the fiducial model, which are partly deviating from the ones in 308 
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Parameter 


Description 


Fiducial value 


Quiescent star formation 
^KS 

Nk 

^crit 


IN Ul lliOjllZjClLIUli Ui XVCllllil^ Lib L ICLVV 

Power-law index in Kennicutt law 
Critical surface density 


1 fi7 V in^'^ M^. Vr\r~^ 

1.4 

6Mq pc-2 


Burst star formation 

/^crit 


Critical mass ratio for burst activity 


0.1 


SN feedback 

^SN 
^rh 

Xrc — infall 


Normalization of reheating fct 

Velocity scale for ejecting gas 
Time-scale for re-infall of ejected gas 


1.3 
2.0 

120 km s-l 
0.1 


Chemical evolution 

y 


Chemical yield 


1.5 


Black hole growth 

/bH, final 
/BH,crit 


Efficiency of conversion of rest mass to radiation 

Mass of seed black hole 

Scaling factor for mass after merger 

Scaling factor for critical BH mass 


0.1 

100 Mq 

0.8 

0.4 


AGN-driven winds 
^wind 


Coupling factor for AGN-driven winds 


0.5 


Radio-mode feedback 

^radio 
^hcat 


Normalization of 'radio mode' accretion rate 
Coupling efficiency of radio jets with hot gas 


2 X 10-3 
1.0 



Table 2. Overview of the different physical mechanisms for galaxy formation assumed in the different SAM versions and implemented 
in the simulations. 



Model 


Gas cooling 


Star formation 


Metals 


Thermal SN feedack 


SN winds 


AGN feedback 


NF 


Yes 


Yes 


No 


No 


No 


No 


SN 


Yes 


Yes 


No 


Yes 


No 


No 


SNWM 


Yes 


Yes 


Yes 


Yes 


Yes 


No 


FULL 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


SIM 


Yes 


Yes 


No 


Yes 


No 


No 



regime, where the cooling radius is larger than the virial ra- 
dius rcooi > ''vir, the cooling rate is set to the gas accretion 
rate, which is governed by the mass accretion history. Slow 
cooling occurs whenever the cooling radius is smaller than 
the virial radius TcooI < '"vir- Here, the cooling rate is calcu- 
lated according to eq. IH The same cooling function is used 
in the simulations, however, the cooling rate is calculated 
locally based on the density and temperature. 

(ii) Photo-ionization: Photo-ionization heating is con- 
sidered in all four SAM versions. It causes halos below a 
certain filtering mass Mf to have a lower baryon fraction 
than the universal average. The collapsed baryon fraction 
as a function of redshift and halo mass is parameterized by 
the expression: 

fb 

/b,coii(^, M.,.) = [i + o.26AfH^)/MviJ3' 



where ft is the universal baryon fraction. The filtering mass 
is a function of redshift and depends on the reionization 
history of the universe (|Kravtsov et al.l |2004| '1. Note that 
photo-ionization heating has very little effect on galaxies 
with circular velocities larger than about 30-50 km/s, and 
therefore plays a minor role in our study, which mainly fo- 
cusses on larger galaxies. In the simulations a UV heating 
background is implemented instead of a filtering mass. How- 
ever, the filtering mass treatment adopted in the SAM is 
based on the results o f numerical hydrodynamic simulations 
(jKravtsov et al.ll2004l ). so we do not expe ct this to intro- 
duce any significant discrepancy (see also iHambrick et al.l 
I2OO9I and references therein). 

(iii) Star formation; The SAM distinguishes between 
quiescent star formation in isolated disks and merger-driven 
starbursts. The quiescent star formation is based on the the 
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empi rical Schmidt-Kennicutt (SK) relation (|Kennicuttl 19891 . 
Il998l ). The star formation rate density is calculated accord- 
ing to 



^SFR 



(8) 



with ^Konn = 1-67 X 10~ , Nk = 1-4, and Egas is the surface 
density of cold gas in the disk. The normalisatio n uses the 
conver sion factor appropriate for a Chabrier IMF l|Chabrierl 
l2003ah . The gas follows an exponential disk (proportional 
to the scale-length of the stellar disk) and only gas above 
a critical surface density threshold Ecrit (~ QMq /pc^) is 
available for star formation. 

Star formation during starbursts is driven by merger 
events. The star formation rate is assumed to be a func- 
tion of the mass ratio and the combined cold gas content 
of the merging galaxies, the bulge to total stellar com- 
ponent and burst timescale. The starburst efficiency as a 
function of these variables is based on hydrodynamic sim- 
ulations of binary galaxy mergers (see references in SOS). 
While in the SAMs we use a 2D implementation for qui- 
escent star formation following the Schmidt-Kennicutt law 
and a simple recipe for starbursts, in simulations, the star 
formation efficiency (in both quiescent and burst modes) is 
determined by the local 3D cold gas density. The normaliza- 
tion of the Schmidt-Kennicutt law was chosen by requiring 
a smooth, isolated disk to lie on the observed relation (see 
ISpringel fc Hernauistl l2003l). 

(iv) Supernova feedback; Exploding supernovae de- 
posit thermal energy in the interstellar medium, which may 
heat the cold gas, and in certain situations may drive winds 
that unbind the gas from the potential well of the dark 
matter halo. In the SAM, these processes are modeled by 
removing cold gas from the galaxy and depositing it ei- 
ther in the hot gas reservoir, where it can cool again fairly 
quickly, or ejecting it from the halo, where it can fall back 
again on a longer timescale. We will refer to the former 
as "thermal" SN feedback and the latter as "SN-driven 
winds". The SN model includes only the thermal SN FB, 
while the SNWM includes both thermal SN feedback and 
SN-driven winds. While supernova driven winds have been 
impl emented in some numerical hydrodynamic simulations 
(e.g. lOppenheimer fc Dav^ lToO^ , only thermal SN feedback 
is implemented in the simulations used in this study. 

The heating rate of the cold gas is given by 



SN 
mrh — (-0 



200 km/s 



(9) 



where and a^h are free parameters and m, is the star for- 
mation rate. We assume the circular velocity Vdisk to be the 
maximum rotation velocity of the dark matter halo, Knax. 
To estimate reasonable values for the paramete r we fol- 
low the prescription of lKauffmann et al.l (|l993l ). in which it 
is assumed that the energy released from supernovae heats 
the gas to the virial temperature of the halo (correspond- 
ing to a value of arh ~ —2). Using this recipe, the thermal 
energy rate is given as 



Et'. 



' ethc 



Qal '7SN i?SN m, — rhrh V^is 



(10) 



where t]sn is the number of supernovae expected per solar 
mass of stars formed (= 4 x 10"'^ Mg^), i?sN the kinetic 
energy of the ejecta from each supernova (~ 10*^ erg) and 



ethormai the efRcieucy with which supernova energy is de- 
posited in the gas, which is highly uncertain. With equations 
[land [in] eo^ is defined as: 



SN 



^thermal 



risN EsN 
(200km/s)- 



(11) 



We assume a value for the thermal efficiency of ethcrmai ~ 

0. 16. 

Additionally we assume to have a kinetic SN feedback, 

1. e. reheated gas can be blown out of the halo. Thereby, the 
fraction of reheated gas, which is ejected from the halo into 
the intergalactic medium (IGM), is given by 



./eject ( Vvir 



14- 



(12) 



with acjcct = 6 and Kjoct a free parameter (« 100 — 
150 km/s). Note that using this definition the total amount 
of released energy from SN explosions is not exceeded. 

Moreover, the ejected gas can re-collapse onto the halo 
at later times and then is available for coo li ng. A s in 
ISpringel et all (|200ll ) and |De Lucia fc BlaizotI (120071 ) the 
rate of reinfall of rejected gas is given by 



'^rcinfall Xrcinfall 



^72cjcet 



dyn 



(13) 



Here, Xreinfaii is a free parameter, rriojoct is the mass of the 



ejected gas outside of the halo and tdy 



rvir/Vvir is the 



dynamical time of the halo. This treatment is quite similar to 
that used in simulations which explicitly include large-scale 
winds (e.g. Oppenhcimcr & Dave 2008), however, as noted 
above the simulations considered here do not include such a 
treatment of winds, and the thermal energy deposition that 
is included is not able to drive gas out of the galaxies. 

(v) Metal enrichment: To track the production of met- 
als, we assume that, together with a parcel of new stars dm* 
a certain mass of metals dMz ~ ydm, is created and instan- 
taneously mixed with the cold gas in the disc. The yield is 
assumed to be constant and is treated as free parameter. 
Whenever new stars are formed, they are assumed to have 
the metallicity of the cold gas at this time step. When metals 
get ejected from the disc due to SN-winds, either the metals 
are mixed with the hot gas or ejected from the halo into 
the 'difi^use' IGM in the same proportion as reheated cold 
gas. Note that only metal enrichment due to Supernovae 
Typell is tracked. Note that in our simulations we consider 
only primordial metallicity cooling and no metal evolution 
is included. 

(vi) Black hole growth and AGN feedback: Every 
"top level" halo (halo with no progenitors) in the merger tree 
is seeded with a black hole with mass ~ IOOMq . Black holes 
can grow by two channels: quasar mode and radio mode. 
The quasar mode is the bright mode of black hole growth 
observed as optical or X-ray bright AGN radiating at a sig- 
nificant fraction o f their Eddington limit (L « (0.1 — l)Z/Edd; 
IVestergaardll2004l : iKollmeier et al.ll2o"oi ). Such bright AGN 
are believed to be fed by optically thick, geometrically 
thin accretion disks ijShakura fc Sunvaevlll97^ ). In contrast, 
AGN activity in the radio- mode is much less dramatic. 
A large fract ion of massive g alaxies are detected at radio 
wavelengths ([Best et al ] |2007l l without showing characteris- 
tic emission lines of clas sical optical or X-ray bright quasars 
jKauffmann et al1l2008l ). Their accretion rates are believed 
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to be a small fraction of the Eddington rate and they are 
radiatively extremely inefRcient. Even if AGN spend most of 
their time in the radio-mode, they gain most of their mass 
during the short and Eddington limited episodes of quasar 
phases which in the model are assumed to be triggered by 
merger events. The energy released during the rapid growth 
of the black holes can drive powerful galactic scale winds 
that sweep cold gas out of the galaxy. 

In contrast to the quasar mode, the radio mode has low- 
Eddington ratio accretion rates, is radiatively inefficient and 
associated with efficient production of radio jets that can 
heat gas in a quasi-hydrostatic hot halo. Assuming Bondi- 
Hoyle ac cretion combined with an isothermal cooling flow 
solution (jNulsen fc Fabianlbood ') we calculate the accretion 
rate in the radio mode 



^^radio — ^radio 



kT 



A(r, Zh) 



M. 



(lO«Me: 



(14) 



Note that the central black hole accretes at this rate when- 
ever hot halo gas is present ('hot mode' accretion, rcooi < 
rvir). The energy that effectively couples to and heats the 
hot gas is given by Lheat = Kheatr/radiriradioC^. Assuming that 
all the hot gas is at the virial temperature of the halo, the 
heating rate is given by 

-^'heat 



W-hcat ~ 



3/4 K 



2 

vir 



(15) 



The net cooling rate is then the usual cooling rate minus 
the heating rate from the radio-mode. Note that in our sim- 
ulations, black hole growth and AGN feedback is not imple- 
mented. 

Note that in merger trees from N-body simulations it 
may happen that the total mass of two merging halos at 
the beginning of a merger event is larger than the mass of 
the merged object afterwards, as during the merger particles 
can become unbound through tidal forces. Therefore, in the 
SAM we impose an upper limit on the hot halo mass of 



Afb 



Mcolc 



Afei 



(16) 



Here, Meject is the mass ejected by winds, Mstar.tot and Mcoid 
are the total star and cold gas masses within the merged halo 
and Mbar is the expected baryonic fraction of the halo. In 
this way, we prevent the sum of all baryonic components in 
the halo from exceeding the universal baryon fraction. 



5 REDSHIFT EVOLUTION OF GALAXY 
PROPERTIES 

In this section, we compare the cosmic evolution of the bary- 
onic components of the galaxies and halos from the direct 
cosmological simulations to the results from the SAMs us- 
ing the dark matter merger trees constructed from the re- 
simulations. Here we only consider the evolution of the cen- 
tral galaxy in the main branch of the merger tree (largest 
progenitor halo) . We divided the 48 halos into three bins ac- 
cording to their halo mass at 2; = (every bin contains 16 ha- 
los) with 4.5 X W^^Mq < Mhaio < 4 x 10"Mq (high -mass), 
1.2 X 10^^ M0 < Mhaio < 4.5 X lO^^M© (intermediate- mass), 
and 2.4 x IO^Mq < Mhaio < 1.2 x lO^^M© (low- mass). All 
comparisons in this Section make use of these bins. 

Note also that a resolution study of the evolution of the 



Redshift 
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2.4«10" M„ < Mj.,„ < 1.2.10" 




2 4 6 B 10 
Lookback time (Gyrs) 

Figure 2. Total baryonic mass Af^ar = A/gas -I- A/», as a fraction 
of the cosmic baryon fraction times the halo mass /bar x Afbaloi 
as a function of lookback time for different semi-analytic models 
(red: no feedback - NF; green: thermal Supernova feedback - SN; 
blue: thermal Supernova feedback, SN-driven winds and metal 
cooling - SNWM; purple: 'full' model including feedback from 
black holes - FULL) and for the SPH-re-simulation (black lines). 
Upper panel: Average values for the high mass bin with 4.5 X 
IO^^A/q < Afbalo < 3.3 X IO^^Mq. Middle panel: Average values 
for the intermediate mass bin with 1.2 X 10^^ Mq < A4aio < 4.5 X 
lO^'^AfQ. Lower panel: average values for halo masses between 
2.4 X 1O"M0 < Mhaio < 1-2 X IQi^M©. 



baryonic component in SAMs based on 2 x and 4 x higher 
resolution simulations for a high- and a low mass halo can 
be found in the Appendix. In both cases the results based 
on the simulations with different resolution are consistent. 



5.1 Baryon fraction 

For a first comparison we compute the total baryonic mass 
A/bar = (Mstar + Afcoid + Mhot) withiu the virial radius of 
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the main halo at every redshift for the simulations and the 
SAMs, respectively. In the re-simulations as well as in the 
SAMs, we consider the hot gas mass within the whole halo 
(i.e. within the virial radius), but the stars and the cold gas 
only of the central galaxy. We neglect contributions from 
substructures, diffuse stars and cold gas and satellite stars 
and cold gas. The baryonic mass is compared to the mass of 
available baryons within each halo, defined as /bar x Mhaio, 
where /bar = 0.169 is the cosmic baryon fraction. 

In Fig. [2] we show the average ratio of Mbar/(,fbar x 
A^haio), as a function of redshift for the three mass bins. The 
simulations are compared to the four SAM variants: NF, SN, 
SNWM, and FULL (for details see section^. We expect the 
simulations to be most directly comparable to the NF or SN 
SAMs, as these SAMs include the same complement of phys- 
ical processes as the simulations. Considering first the NF 
model, we see that at low redshift, the SAM overestimates 
the baryon fraction in high mass halos, nearly agrees in in- 
termediate mass halos, and slightly underestimates it in low 
mass halos. At high redshift, the NF SAM overestimates 
the baryon fraction at high redshift in high and intermedi- 
ate mass halos, by a somewhat larger factor in the former. 
Turning next to the SN SAM, we see that the SAM pre- 
dicts baryon fractions that are everywhere higher than the 
simulation results, though much more so for the high and 
intermediate mass halos. This is because the "thermal" SN 
feedback removes baryons from satellite galaxies, which are 
not counted in this census, and deposits them in the hot 
gas component which is included here. In the SNWM and 
FULL model, we see the impact of the SN-driven winds, 
which remove baryons from low-mass halos. The additional 
metal cooling in the SNWM and the FULL model does not 
have any impact on the total baryon fraction (only on the 
individual components, see next sections.) AGN feedback 
(in the FULL model) mainly prevents hot gas from cooling, 
so does not affect the total baryon fraction significantly, but 
will be important for the fraction of stars and hot gas, which 
will be discussed later. 
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5.2 Cold gas and stars 

In Fig. 13] we show the evolution of the mass of condensed 
baryons (stars and cold gas, M* -|-Mcoid) as a fraction of the 
total baryon mass as a function of redshift for the three mass 
bins. There is fairly good agreement between the simulations 
and the NF model over the whole redshift range, although 
the SAM is a little low at high redshift and a bit high at low 
redshift, particularly in the high and intermediate mass bin. 
In the SN model, the condensed baryon fraction is lowered 
by an almost fixed factor relative to the NF model, and is 
significantly lower than the simulation results. This suggests 
that the "thermal feedback" implemented in the simulation 
is less effective than that included in the SAM. 

It may seem curious that the SNWM model results are 
higher than the SN model, in fact close to the NF model 
in the high and intermediate mass bin. This is because the 
SNWM model includes metal cooling, leading to more ef- 
ficient star formation, while the SN model does not. The 
enhanced cooling rates partly compensate for the removal 
of cold gas via the SN-driven winds. In particular for the 
high- and intermediate mass bin the effect of the SN winds 
is almost the opposite of that of metal cooling, resulting in 



Figure 3. The evolution of the mass of the condensed baryons 
M, + Mcoid as a fraction of the total baryon mass as a function of 
redshift for the three mass bins. The condensed baryon fraction in 
halos of all masses is in reasonable agreement with the NF model 
at all redshifts. We can see that SN FB most strongly affects the 
low-mass bin, while AGN FB affects the high mass bin. 



condensed baryon fractions close to the NF model. Only in 
the low-mass bin, the effect of SN winds is not completely 
compensated by metal cooling (so condensed baryon frac- 
tions are below the NF model). 

Finally, we see that in the FULL model, the AGN FB 
begins to quench star formation in the massive halos after 
about 2 ~ 2, while it has little effect on the lower mass bins. 

In Fig. |4] we plot the evolution of the mean cold gas 
fraction of the central galaxy (in Fig. [3] we plotted cold gas 
plus stars). The efficiency of the conversion of cold gas to 
stars is clearly very different between the simulations and 
SAMs. In both cases (simulations and SAMs), the final cold 
gas fraction is increasing with decreasing halo mass. For 
the NF model, the SN and the SNWM model the cold gas 
fraction varies only slightly over time and is significantly 
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Figure 4. Evolution of mean cold gas fraction Mcoid/(/bar x 
Mjiaio) of central galaxies. Mass bins and colors are the same as in 
Fig. [2] In all simulations the cold gas is depleted more efficiently 
than in the SAMs due to tlie large star formation efficiency at 
high redshifts. The red dotted line shows the cold gas fraction 
assuming ten times higher efficiency for star formation in the NF 
model (see Eg. 1171 1. 



higher (about an order of magnitude since 2 = 2) than for 
the simulations. This shows that the inclusion of SN FB 
has little impact on the gas fractions of galaxies. Only the 
FULL model shows a much stronger decrease of the gas frac- 
tion with cosmic time for massive galaxies due to the radio 
mode feedback. The initial cold gas fraction, at high red- 
shifts 4 < 2 < 8, is almost the same for the simulations 
and SAMs. With evolving cosmic time the cold gas content 
decreases more rapidly in the simulations due to the more ef- 
ficient conversion into stars. The cold gas in the simulations 
is already converted into stars at high redshift and there is 
almost no more cold gas left to turn into stars at lower red- 
sh ifts. This is simi l ar to the results found in the comparison 
of ICattaneo et al.l l|2007l '). 



Figure 5. Comparison of the stellar baryon fraction of the central 
galaxy between simulations and SAMs. The mass bins and colors 
are the same as in Fig. [2] For all mass bins the simulations agree 
best with the NF SAM, but form significantly more stars at high 
redshifts z > 1. We can reproduce the simulation results fairly 
well with the NF SAM if we increase the star formation efficiency 
parameter by a factor of ten (red dotted lines). At high masses 
the AGN feedback (FULL) and at low masses the SN feedback 
(SN and SNWM) reduce the stellar baryon fractions in the SAMs. 



In Fig.[5]we show the corresponding fraction of available 
baryons that are converted into stars in the central galaxy 
M,/(fha. r X Mhalo), sometimes termed baryo n conversion ef- 
ficiency (|Guo et al.ll2009l : lMoster et al.ll2010l '). In general, all 
simulations predict a decreasing (high-mass) or almost con- 
stant conversion efficiency with redshift (low-mass), whereas 
most SAMs predict increasing conversion efficiencies with 
the exception of high- mass galaxies in the FULL model with 
AGN feedback. 

At low redshift z < 0.6 the conversion efficiencies agree 
well between the simulations and the NF model, with higher 
values for lower mass galaxies. However, at high redshifts 
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z > 1 the conversion efficiencies are significantly liiglier 
for tfie simulat i ons. T fiis is in contrast to the results of 
ICattaneo et al.l ()2007[ l. where the stellar masses agree at 
high redshift, but the SAM masses are larger than in sim- 
ulations at low redshift. The difference in the behaviour of 
the SAMs and the simulations can be explained in terms of 
star formation efficiency. We changed the normalization of 
the SK-relation in the SAM by introducing a factor r, in 
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AkS ^Nk 



(17) 



with r* « 0.1. The results of the NF model with this elevated 
star formation efficiency for the stellar and cold gas mass 
evolution is shown in Figs. [4] and [5] At high redshift, gas is 
more efficiently depleted and converted into stars, resulting 
in a better agreement between this 'high SFE' model and the 
simulations. However, for the high-mass and intermediate- 
mass bin, the 'high SFE' model overpredicts the stellar frac- 
tion a t low redshifts, suggesting that the ICattaneo et al.l 
|2003) SAMs may also have had a higher SFE, and this could 
explain the discrepancy between their results and our initial 
results. In all three mass bins, the NF model produces the 
most massive stellar components. Again, the more efficient 
cooling due to metals in the SNWM and the FULL model 
is cancelled by the effect of winds and, for massive galaxies, 
also by AGN feedback, resulting in lower stellar masses than 
in the NF model. The separate behaviour of metal cooling 
and SN winds is the same as for the total condensed baryon 
fraction. 



5.2.1 Star Formation Rates 

To confirm the previous findings we compare the star forma- 
tion rates in Fig. [B] At very high redshifts z > 4, the SFRs 
in the simulations are much higher than in the SAMs. Only 
by assuming more efficient star formation in the NF SAM 
(r* — 0.1) do we obtain a reasonable match to the simula- 
tions. However, at z < 1.5, the high SFE model results in 
similar SFRs as the original NF model, as the larger SF effi- 
ciencies at high redshift lead to a more rapid depletion of the 
cold gas. In the simulation, the cold gas is rapidly turned into 
stars, resulting in lower SFRs at low redshifts compared to 
SAMs because of gas depletion. Only the FULL model shows 
a strongly decreasing SFR with decreasing redshift due to 
radio mode feedback, which becomes especially important 
for low redshifts and lar ge halo mas s es. T his result is con- 
sistent with the study of Saro et al] (l2010|), who compared 
their stripped-down versions of SAMs (with no feedback) 
to simulations, and found find higher SFRs in the simula- 
tions for all galaxies within a cluster (central and satellites) 
at hig h redshifts and lower SFRs at low redshifts. In ad- 
dition, IStringer et al.l |2010) find a similar discrepancy for 
the specific star formation rates at high redshifts (larger in 
simulations than in their SAM) and good agreement for low 
redshifts. 

To better understand this discrepancy between sim- 
ulations and SAMs we take a closer look at the re- 
spective implementat i ons of star formation. According to 
ISpringel fc Hernguistl (|2003l ). stars in the simulations are 
formed locally out of cold gas with the star formation rate 
density proportional to the local three-dimensional density 
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Figure 6. Evolution of the star formation rates in simulations 
and SAMs for three mass bins. The mass bins and colors are 
the same as in Fig. [2] At high redshifts, SFRs are higher in the 
simulations than in the SAMs, leading to more rapid depletion of 
the cold gas and very low SFRs at low redshifts. 



of gas to the power of 1.4, psF oc p^'^/t,. The star for- 
mation timescale t« was set to reproduce the observed lo- 
cal Schmidt-Kennicutt relation (SK) for a simulation of a 
smooth, isolated disk-dominated galaxy set up to resemble 
the Milky Way. In the SAMs the cold gas is assumed to settle 
into smooth exponential disks and stars form according to 
the SK-relation, implemented in terms of surface densities. 

In Fig. [3 we plot the SFR surface density versus the 
surface density of the cold gas for the simulated galaxies at 
z — 2 and z = 4 within 1/10 r^ir for all re-simulations. The 
black da shed l i ne is the SK-relation assuming a Salpeter 
IMF (Salpeted [l955h . as given in the original Kennicutt 
papers and as implemente d in the simulations following 
ISpringel fc Hernguistl (|2003D . We would naively expect the 
simulations to follow this line. Th e red solid line sh ows the 
SK-relation for a Chabrier IMF (|Chabrieij [2003bl ), as as- 
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Figure 7. Star formation rate surface densities versus cold gas 
surface densities for the simulated galaxies within 1/lOrvir- Black 
and blue stars correspond to different re-simulations at 2 = 2 or 
2 = 4, respectively. The red solid line illustrates the Kennicutt 
relation implemented in the SAM with a normalization consistent 
with a Chabrier IMF, while the black dashed line the o ne used 
to normalize the simulations llSpringel fc Hernauistll2003l '). which 
assumes a Salpeter IMF. 



sumed in the SAMs. At a given gas surface density, the 
SFR surface densities of the simulations he mostly above 
the expected SK-relation. The change of normahzation as- 
sociated with converting from Salpeter to Chabrier cannot 
account for the increased star formation efficiency in the sim- 
ulations. In general, star formation in the cosmological sim- 
ulations is about a factor of five more efficient than for sim- 
ulat ions of smooth isolated disk s using the identical model 
fsee lSpringel fc Hernauistll2003f ). This discrepancy is a con- 
sequence of the clumpy structure of cold gas in the cosmolog- 
ical simulations. In the clumps the gas can reach higher local 
densities than in the idealized s moot h disks that have been 
used by ISpringel fc Hernguistl (|2003l ') to calibrate the star 
formation timescale by matching the SK-relation. As the im- 
plemented SK-relation is not linear, the structure of the cold 
gas distribution plays an important role for t he overall star 
forma tion efficiency within the galaxies (see iTevssier et al.l 
l|2010l ) for a discussion on galaxy mergers). In other words, 
for any star formation model with a non-linear dependence 
on the local gas density (exponent larger than unity), a more 
clumpy gas distribution will effectively increase the star for- 
mation efficiency. These combined effects explain the much 
higher SF efficiencies at high redshift in the simulations rel- 
ative to the SAMs. 



5.2.2 Modes of Stellar Mass Growth 

In the hierarchical picture, galaxies can grow their stellar 
masses in two ways: 1) by converting cold gas into stars 
in situ 2) by accreting already formed stars via mergers. 
We refer to these two modes as "in situ" and "accreted". 
Our simulations exhibit two phases of growth, with a rapid 
early phase at z > 2 during which stars are formed in situ 
from infalling cold gas, followed by an extended phase at 
z < 3 during which the growth is pri marily due to accr etion 
of stars formed in external galaxies (|Oser et al.ll201ol ). We 
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Figure 8. Fraction of in situ to accreted stellar mass versus red- 
shift for different halo masses. The mass bins and colors are the 
same as in Fig. [2] In situ star formation dominates over accretion 
in all the SAM variants at all redshifts, in contrast to the simu- 
lations for which accretion dominates at late times, especially in 
massive halos. 



now investigate whether the SAMs show the same behavior. 
In Fig. |8l we show the fraction of cumulative in situ over 
accreted stellar mass as a function of redshift for the three 
different mass bins. For the SAMs the qualitative trend of 
a decreasing fraction of in situ growth is reproduced for the 
high mass bin. However, the fraction of in situ formed stars 
dominates over accreted stars for all models, all masses, and 
at all redshifts. This is in contrast with the simulations, 
where accretion dominates over in situ f ormation for massi ve 
systems at low redshifts as discussed in lOser et all (|2010h . 

We note several interesting trends in the in situ to ac- 
creted fraction as we vary the physics in the SAMs. Adding 
thermal SN FB increases the in situ fraction at all redshifts 
and in all mass bins. This is presumably because it sup- 
presses star formation in low-mass satellites which are the 
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Figure 9. Comparison of the mean in situ (left column) and accreted (right column) stellar masses in different mass bins. The mass 
bins and colors are the same as in Fig. [21 The in situ stellar masses in SAMs agree with the ones in simulations reasonably well, whereas 
the accreted stellar mass is smaller in the SAMs than in the simulations. 



source of accreted stars. Adding the SN-driven winds and 
metal cooling further increases the in situ fraction, again at 
all redshifts below z ~ 4. Switching on AGN FB increases 
the in situ fraction at high redshift and decreases it at low 
redshift in the high mass bin (and to a lesser extent in the 
intermediate mass bin). This is because the radio mode feed- 
back shuts off cooling at late times in massive halos, remov- 
ing the supply of new gas needed to fuel ongoing in situ 
star formation. Interestingly, increasing the star formation 
efficiency in the NF model has almost no effect on the in 
situ to accreted fraction. This is presumably because the SF 
efficiency is increased in the central and (accreted) satellite 
galaxies alike. However, if the SFE were higher in high red- 
shift galaxies than at low redshift, this would presumably 
increase the accreted fraction in present day galaxies. This 



may be part of the reason for the higher accreted fractions 
in the simulations. 

In Fig. [9l we show the evolution of the cumulative mass 
of insitu and accreted stars separately for the simulations 
and the various SAM variants. Here we can see that the NF 
SAM actually reproduces the growth of m situ stellar mass 
fairly well, though overproducing the in situ mass at low red- 
shift somewhat, especially in the highest mass bin. We spec- 
ulate that gravitational heating in the simulations prevents 
some of the late cooling in the highest mass bin a nd leads 
to lower in s i tu stellar mass than the NF SAM (Naab et al.l 
l2007l . l2009al : [Johansson et al.ll2009l : iFeldmann et al.ll2010l ). 
The radio mode AGN FB in the FULL model leads to a 
similar suppression of this in situ mass growth in the mas- 
sive halos. The NF model with increased SFE gives an even 
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better match to the simulations at high redshift. The dis- 
crepancy arises from the much lower accreted masses in the 
SAM. Here again, the SF model with high SFE comes the 
closest to matching the simulation results, but it still falls 
short by a considerable amount. 

Part of the reason for the lower predicted accreted 
mEisses in the SAMs is that the SAMs used here only al- 
low cooling onto the central galaxy in the halo, effectively 
assuming that the hot gas reservoir of a satellite galaxy 
is stripped as soon as it enters the virial radius of the 
host. This is known to result in satellites that are too red 
and have star for mation rates that are too low compared 
with observations (|Kimm et al.ll2009l ). It will also truncate 
their star formation, resulting in a smaller amount of stel- 
lar mass that will eventual ly be accreted when they merge 
IIKhochfar &: Ostrikejiiooi) . 

5.3 Hot halo gas 

The evolution of the mean hot gas fraction is shown in 
the three panels of Fig. [10] In all but the SN SAMs 
the hot gas fraction increases with increasing halo mass, 
which is qualitatively similar to the simulations. In the 
simulations this effect is caused by shock heating of in- 
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SAMs, except for the SN model, where the supernova energy 
input heats most of the available gas to the virial temper- 
ature of the halos, keeping the hot gas fraction constant 
independent of halo mass. For the SNWM and FULL mod- 
els, the supernova winds drive some of the hot gas out of 
the low mass halos. In contrast, the additional metal cooling 
shows a negligible effect on the evolution of the hot halo gas. 
The additional effect of Radio mode heating (FULL model) , 
which prevents late cooling in massive halos and therefore 
leads to larger amounts of hot gas, is apparent for the in- 
termediate and high mass galaxies. The NF model agrees 
fairly well with the simulations in all mass bins 
but substantially overpredicts the amount of hot gas at high 
redshift. 

To understand the differences in the hot gas content at 
high redshift we investigate the gas accretion modes onto the 
central galaxies. For the simulations we distinguish between 
hot and cold accretion by considering the highest tempera- 
ture a gas particle had before it was accreted onto the galaxy, 
i.e. 1/lOth of the virial radius. We use the same definition 
to distinguish betwee n between hot and c old gas as given by 
equation |3] (similar to lKeres et al.l[2005bl ). In the SAMs, we 
distinguish between hot and cold mode accretion (slow and 
rapid cooling) depending on whether the ratio of the cooling 
radius to the virial radius rcooi/ rvir is larger (cold mo de) or 
smaller (hot mode) than unity (jWhite fc Frenklll99lh . The 
distinction between hot and cold mode accretion approxi- 
mately specifies whether the gas was heated to the virial 
temperature of the host halo before it was accreted onto the 
galaxy (hot mode) or was directly accreted without being 
heated (cold mode). The cold mode accretion is meant to 
repre sent the 'cold flows' recently discussed in the litera- 
ture l|Keres et al.ll2005bl . l2009l : lOser et aLllioiol : iDekel et all 
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Figure 10. Evolution of the mass fraction of hot gas in simu- 
lations and SAMs normalized to the available mass in baryons. 
The mass bins and colors are the same as in Fig. [2] The SAMs 
tend to overpredict the mass of gas at high redshift relative to 
the simulations. 



l2009h . Note that for the SNWM and FULL SAMs, we have 
not substracted the heating rates or rates of blown-out gas 
due to feedback processes from the accretion rates. 

In Fig. [TT] we show the comparison of the total (left 
panels), hot (middle panels) and cold mode (right panels) 
gas accretion rates onto the central galaxies as a function 
of redshift. For all SAMs, the total gas accretion rates onto 
the galaxies are significantly higher than for the simulations. 
This is caused solely by higher hot mode accretion rates from 
the generally larger hot gas reservoir in particular at high 
redshift (middel panel in Fig. 1111 and see Fig. llOp . The cold 
mode accretion rates are much lower in the SAMs than in the 
simulations, and in the SAMs without metal cooling (which 
are more relevant to compare with these simulations), the 
cold mode is truncated at z 4 for massive halos, z ^ 1-1.5 
for low mass halos, while it declines smoothly until 2 in 
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Figure 11. Comparison of the mean accretion rates of all gas (left column), the rate of hot mode accretion (intermediate column) and 
of cold mode accretion (right column) in simulations and SAMs onto the central galaxies. The mass bins and colors are the same as in 
Fig. [2] Note that in the FULL and SNWM we have not substracted heating rates or rates from blown-out gas from the accretion rates. 
Compared to simulations, hot mode accretion is overestimated in all SAMs, whereas cold mode accretion is in general underestimated. 



the simulations. However, the predicted rates of cold mode 
accretion are much higher in the SAMs that include metal 
cooling. 

In contrast to our results, ICattaneo et"ai] (|2007l 1 find 
a reasonably good match for the evolution of the hot gas 
content as well as for the hot and cold mode accretion rates 
in their SAM version without any feedback. Ifowever, they 
inclu de metal cooling in their SAM, but not in their simula- 
tions. iBenson fc Boweil (|201ll ') compare cold and hot mode 
accretion rates from SAMs to simulations, varying the super- 
nova feedb ack and conditions for the r apid cooling regime ac- 
cording to lBirnboim fc Dekej (|2003al ). They concluded that 
cold-mode physics is already adequately accounted for in 
SAMs — but they also used simulations with only primor- 
dial H fc He cooling, but included metal cooling in their 
SAMs. lLuet"ZI(|201(]| ) assume only H fc He cooling in their 
SAMs OS well as in the simulations (like we do) and find 
qualitatively similar results to ours: a discrepancy for the 



hot halo gas fraction at high redshift associated with larger 
hot mode and smaller cold mode accretion rates in the SAMs 
than in the simulations. 

This sugg e sts t hat the agreement presen ted in 
ICattaneo eralT l|2007f ) and iBenson fc Boweij (120111 ) is for- 
tuitous, and arises because of the enhanced cold flows re- 
sulting from the metal cooling included in the SAMs. It is 
unclear whether this agreement would persist if metal cool- 
ing were also included in the simulations, but it appears that 
the agreement is not nearly as good as they claim when 
metal cooling is omitted from both techniques. Therefore, 
provided that the physical mechanisms taken into account 
in SAMs and simulations are the same, we might conclude 
that the difference between cold and hot mode accretion 
rates in SAMs and simulations is a general result, indepen- 
dent of the specific SAM that is considered. This indicates 
that the cooling recipe should be improved in SAMs. In ad- 
dition, the recipes for gas accretion in SAMs do not currently 



Semi- analytic models and cosmological simulations 17 




10.5 11 11.5 12 12.5 13 13.5 10.5 11 11.5 12 12.5 13 13.5 10.5 11 11.5 12 12.5 13 13.5 14 
log(M,,,,/Me) log(M,,,yMe) log(M,„yMe) 



Figure 12. Stellar mass versus dark matter halo mass for simulations (upper row), NF model (middle row) and the FULL model (lower 
row). The left column shows the dependence on halo mass for z = 0, the middle one for z = 1 and the r i ght on e for z = 2. The black lines 
with the light grey shaded areas sh ow the fit for the halo occupation distribution from lMoster et al.l | |2010| '|. The green lines with the 
dark grey shaded areas the one from lWake et al. Il l201ll 'l. The difi^erent black, dashed and dashed-dotted lines illustrate the total expected 
baryon fraction (/bar * ^^halo): 50% and 10% of the total expecte d ba ryon fraction. The black and green thin, vertical lines show the 
observational mass limit of the survey used in lMoster et alj l l20ld) and lWake et al.l 1I2OIII '). respectively. 



allow co-existing cold_and hot gas accretion as seen in simu- 
lations. For this. lLu et a l. (2010) proposed a new model that 
explicitly incorporates cold-mode accretion independent of 
the hot halo gas. By fitting the hot and cold gas fraction 
in simulations as a function of redshift and halo mass and 
assuming accretion onto the galaxy within a free-fall time 
they calculate the accretion rate of the cold component and 
thus, achieve a better match of their SAMs to simulations. 



6 COMPARISON TO OBSERVATIONS 

In this section we compare our results from the SAMs and 
simulations to observational data and empirical constraints 
at different redshifts. We focus on two key observational 
constraints: the relationship between halo mass and stellar 
mass (the Mgai — Mhaio-relation) and the relationship be- 
tween stellar mass and star formation rates (rhstar — Mgai- 
relation) . 



Fig. [12] shows the relation of galaxy mass and dark 
matter halo mass for z — (left panels), z — 1 (mid- 
dle panels), and z — 2 (right panels). We show the NF 
and FULL SAMs, and the simulations. We also show the 
empirical constrain ts on the A^gai — A/haio-relation from 
iMoster et al] (|2010l ). which were obtained by asking how 
halos and sub-halos in a N-body simulation must be pop- 
ulated in order to reproduce the observed stellar mass 
functions at different redshifts (halo abundance match- 
ing). The thin, black vertical lines illustrate the observa- 
tional lower observation al mass limi t . Also shown are the 
similar constraints from IWake et al.l l|201lh . which are de- 
rived using galaxy clustering data from the NEWFIRM 

Medium Band Survey between 1 < z < 2 (see also 

Wechsler et"an,2 006: Z heng et al.ll2007l: IConrov fc Wechslei 



2OO9I : IGuo eTai.! 120091 : IZehavi et al.l l20l'o[ 
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201c ). Note that the fit ting functions of iMoster et all (I2OI0I ) 



and Twake erahl l|201ll) are somewhat different, in particu- 
lar at the low mass end (lower observational mass limit is 
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shown by thin, green hnes). This is not surprising, as they 
were derived using different methods and from different ob- 
servational data sets. As our sample consists of only 48 ha- 
los covering a mass range between approximately 10^^ Mq 
and 10^^ Mq, the comparison to observational data is not 
statistically rigorous due to the relatively small number 
of simulated halos. However, the 48 halos have been cho- 
sen randomly from a cosmological simulation and therefore, 
they should generally follow the mean abundance matching 
trends that relate halo mass to stellar mass. At all redshifts, 
the simulations overpredict the stellar masses at a given halo 
mass by about a factor of two for halos more massive than 
IO^^Mq. At higher redshifts the progenitor galaxies have 
lower masses, and deviate more from the expected distri- 
bution. Kt z = 2 the difference can be almost two orders 
of magnitude for halos of ~ 10^^ M©, in line with the find- 
ings of the previous section — at high redshift, gas is very 
efficiently converted into stars in the simulations. Implemen- 
tation of more efficient feedback from supernovae would help 
to solve this problem. Indeed it has been shown that simula- 
tions that do include effective SN FB a gree much bettor with 
expectations (see e.g. Scannapieco et al. 2009: Sawala e t al.l 
l2010l : lGenel et al.ll2010l : lGovernato et al.ll201ol and references 
therein). In addition, it is apparent the simulations require 
an additional process that can quench star formation at late 
times in massive halos, such as radio mode AGN feedback. 

For high mass halos, the NF model (middle row) pre- 
dicts galaxy masses close to the relation &t z = 1 and z = 2 
due to less efficient star formation at high redshifts (see Fig. 
[SJ, but, like the simulations, overpredicts the stellar masses 
in low mass halos. By z = the offset is about as large 
as for the simulations. For the FULL SAM there is good 
agreement at 2 = 0, which is not surprising because the 
model was tuned to match the observed stellar mass func- 
tion. However, the FULL model still predicts galaxy masses 
that are about a factor of two to three too high for low mass 
halos (log(Mhaio) < 11.5) at high redshift. This is related to 
the excess of low mass galaxies at high redshif t and other 
connected problems discussed in lFontanot et al] (|2009l ). and 
is seen in both SAMs and hydro simulations from several 
groups. It is likely that these problems are due to limita- 
tions in our current understanding or implementation of the 
physics of star formation and/or SN feedback. 

The relation between the star formation rates and the 
galaxy stellar masses Mstar is shown in Fig. [TH] for redshifts 
z = 0, 2 = 1, and 2 = 2. We distinguish between star- 
forming and non-star-forming galaxies (illustrated as crosses 
or open squares, res pectively) using a criterion according 
to lFranx et al.l l|2008l ): galaxies with specific star formation 
rates SFR/Mst!a smaller than 0.3 x t^ubbiei ^'^^ considered 
to be quiescent, whereas galaxies with larger specific star 
formation rates are assumed to be star-forming. The black, 
solid lines always refer to the obs erved relation at th e corre- 
sponding redshift (2 = 0: SP SS, lElbaz et all (|2007l); 2 = 1: 



GOOD S, lETbaz et all l|2007l ): 2 = 2: GOODS. iDaddi et all 
l|2007h ) f or star-forming gala xies. Note that we show the re- 
lation of iDaddi et all l|2007h at 2 = 2 re-normalised 0.3 dex 
downwards, following the re-calibration of SFR derived from 
24 micron luminosit y based on recent Herschel observations 
l|Nordon et al.ll2010l ). We show Fig. [T3]in order to illustrate 
the variation in the population of the SFR- Mstar plane for 
different implementations of the SAMs and the simulations. 
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Figure 13. Star formation rate versus stellar mass for simulations 
and two different semi-analytic models (NF and FULL). Solid 
line s illustrate the obs erved re lations at different redshifts: 2 = 
0, l: lElbaz et"all l|2007l ). 2 = 2: IPaddi et all l|2007l V Dashed lines 
illustrate the corresponding 1 — cr-range of the scatter. Note that 
in light of the recalibration of SFR from 24 /^m observations based 
on Herschel results (see text), we have plotted the Daddi et al. 
relation shifted downwards by 0.3 dex. 



Note that this is not intended to be a quantitative com- 
parison to observations, but only a qualitative illustration. 
In particular, we can see that — with regard to the SAMs 
— only in the FULL model, do we have SF and non-SF 
(quenched) galaxies co-existing as in the observations. This 
reffects the well-known need for some mechanism, such as 
AGN feedback, to quench SF in massive galaxies. 

In general, the simulations underpredict the star for- 
mation rates for star forming galaxies at all redshifts, but 
most notably at 2 > 1, due to the high star formation ef- 
ficiencies at even higher redshifts and the resulting gas de- 
pletion (see section [5^ . Once again, implementing more ef- 
fective supernova feedback would presumably suppress star 
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formation in the small, high redshift progenitors of these 
gal axies and result in higher SFR at these reds hifts (see 
e.g. lQppenheimer fc Dav ^20081 : iGenel et al.ll2010l ). The NF 
SAMs fit the SF sequence at z = 0, but have too many high 
meiss galaxies with high SFR. In the FULL model, these 
galaxies are quenched by r adio mode feedback, in agree- 
ment with observations (see ISomerville et al.ll2008bl ). Both 
the NF and FULL SAMs SF sequences are about a factor of 
two too low at 2 ~ 1 and 2. This also seems to be generic to 
many SAMs, as shown by iFontanot et alj l|2009l ) and oth- 
ers. It is interesting to note that star forming galaxies in 
both the NF and FULL SAM galaxies all lie on the same 
SF sequence. This implies that SN feedback simply moves 
galaxies along the sequence, reducing both the stellar mass 
and the SFR such that the galaxies remain on the same rela- 
tion, and is due to the self-regulating nature of SN feedback 
in the SAMs. In constrast, AGN FB quenches star formation 
and moves massive galaxies off of the SF sequence (seen in 
the FULL model). 



7 DISCUSSION AND CONCLUSIONS 

In this paper we present a detailed comparison between a set 
of 48 cosmological hydrodynamic zoom simulations and dif- 
ferent stripped-down versions of semi-analytic models based 
on dark matter merger trees extracted from the simulations. 
The hydro simulations are run using the entropy-conserving 
formulation of SPH with the Gadget-2 code, and include 
atomic cooling assuming a primordial composition of H and 
He, a UV background radiation field, star formation, and 
supernova feedback. We present results from a "no feed- 
back" (NF) version of the SAM which contains the same 
physical ingredients as the simulations, and also versions 
that include thermal feedback by SNae (SN), a version that 
also includes large-scale SN-driven winds and metal cool- 
ing (SNWM), and the FULL version which includes all the 
previously mentioned ingredients as well as AGN feedback. 
With this approach we can compare the predictions of the 
two methods over two orders of magnitude in halo mass at 
unprecedented resolution. 

The two approaches try to answer the same questions 
but differ in methodology. In the simulations the full dy- 
namical and hydrodynamical evolution of the systems is fol- 
lowed by solving the equations of motion computationally. 
Additional physical processes like star formation and super- 
nova feedback are included using sub-resolution models. The 
semi-analytic models are based on the computed dark mat- 
ter accretion history and approximate the gas physics, star 
formation, and feedback processes with simplified recipes. 
Our analysis is focused on the cosmic evolution of the baryon 
content in the central galaxies of the main branch of the 
merger trees and its division into various components (stars, 
cold gas, and hot gas), as well as how those galaxies acquired 
their gas — whether through "cold" or "hot" mode accretion 
— and their stellar mass (e.g. through in situ star formation 
vs. accretion). 

The results of our comparison are quite rich, with some 
surprising agreement and some striking disagreement. First, 
we note that we expect the results of the simulations to 
lie somewhere in between the NF and the SN SAMs, since 
these SAMs include the same physical processes as the sim- 



ulations. In most cases, the agreement is best between the 
simulations and the NF model, suggesting that the SN FB 
implemented in the simulations has little effect. The NF 
SAMs produce very good agreement with the simulations 
for the mass of cold gas plus stars at all redshifts and for 
all halo masses. The SAMs slightly underestimate this "con- 
densed baryon" fraction at high redshift (z > 1) and over- 
estimate it at low redshift. This indicates that the overall 
cooling and accretion rates in the SAM and the simulations 
must be similar. The NF SAM also produces fairly good 
agreement (better than ~ 20% since z ~ 4) with the overall 
baryon fractions (i.e. hot gas plus cold gas and stars in the 
central galaxy) in the simulations, here overestimating the 
baryon fractions at high redshift in high and intermediate 
mass halos. 

A striking difference is that when we study the evolu- 
tion of the stellar and cold gas components separately, we 
find that the cold gas fractions agree at very high redshifts, 
but the gas is consumed much more rapidly in the simula- 
tions, leading to cold gas fractions at all redshifts less than 
about z ~ 3-4, and in halos of all masses, that are much 
lower (by up to two orders of magnitude) than in the NF 
SAM. Correspondingly, we find much higher stellar masses 
in the simulations than the SAMs at high redshift, although 
they converge to almost the same value as the NF SAMs at 
z = 0. We interpret this as an indication that the star for- 
mation efficiency is much higher in the simulations than in 
the SAMs, and the reason for the convergence in the stellar 
masses at low redshifts is because nearly all available gas 
has been consumed in the simulations. This conjecture is 
supported by our finding that if we boost the star forma- 
tion efficiency in the NF SAM by a factor of ten, we find 
excellent agreement with the stellar mass fraction evolution 
in the simulations, and improved agreement with the cold 
gas fraction evolution. 

However, both the simulations and SAMs supposedly 
adopt the same empirical Schmidt-Kennicutt relation be- 
tween cold gas density and star formation rate. How can 
the star formation efficiencies be so different? We note sev- 
eral differences in the implementation of the star formation 
recipe in the simulations and SAMs. In SAMs, the only avail- 
able information about the structure of the star forming gas 
in galactic disks is an estimate of the scale radius of the 
total baryonic component of the disk, which comes from an- 
gular momenturn coii servation arguments (|Mo et al.lll998l : 
ISomerviUe et aLlbOOSl ). The SAMs then make a series of as- 
sumptions — that the gas is in a smooth, thin exponential 
disk with a radius that is a simple multiple of the stellar 
scale radius — and apply the Kennicutt relation in terms 
of the predicted gas surface density. Only gas above a criti- 
cal surface density is allowed to form stars. In constrast, the 
simulations provide detailed 3D predictions for the structure 
of the cold gas in galaxies, and implement the SK relation in 
terms of 3D volume density (also applying a threshold for SF 
in terms of a critical volume density). It is well known that 
high redshift galaxy assembly in cosmological simulations 
is dominated by clumpy, high density, cold mode accretion. 
Disks may be more compact than in the idealized case of 
perfect conservation of angular momentum, and are thick 
and clumpy. The adopted SF recipe is super-linear in the 
gas volume density (i.e. the exponent in the SK relation is 
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larger than unity), and therefore star formation will be more 
efficient in a clumpy gas distribution than in a smooth one. 

The appropriate values of the SF and SN feedback effi- 
ciency parameters for the simulations were obtained by tun- 
ing them to match the observed Kennicutt relation for an 
idealized, smooth thin exponential di sk, designed to resem- 
ble a Milky Way-like galaxy aX z — Q (ISpring el fc Hernguistl 
l2003f ). Using these same values for the parameters, we find 
that the high redshift galaxies in our cosmological simula- 
tions lie about a factor of five above t he Kennicutt relation 
with the normalization adopted by ISpringel fc Hernguistl 
(|2003l ). A second, more minor issue, is that the SAMs 
were tuned to match a Kennicutt relation with a nor- 
malization a factor o f two lower than the one used by 
ISpringel fc Hernguistl (j2003l l. reflecting a different choice of 
IMF in the conversion from observed flux to SFR. 

This seems to explain the reason for the factor of ^ 10 
higher SFE in the simulations relative to the SAMs, but 
begs the question: is this a robust prediction of the sim- 
ulations that should be taken seriously? Are these higher 
star formation efflciencies in high redshift galaxies really 
physical? There are several issues that are relevant here. 
First, the predicted "clumpyness" of the disks is highly sen- 
sitive to the assumed sub-grid recipes. For example, im- 
plementation of more effective SN FB would reduce the 
clumpyness of the disks at high reds hift, but might in - 
crease the clumpyness at z « 2(e.g. iGenel et al.l I2OI0I ) 
in the simulations. Observed disks at high redshift are 
known to be more clumpy than nearby o nes ("Gen zel et al.l 
I2OO6I . l2008l . I2OIOI : iForster Schreiber et al.l 2009. 20l"ll). but 
it remains highly uncertain which recipe for SN FB will 
produce the "correct" degree of clumpyness and overall 
structure for statistical samples of high redshift galaxies 
while simultaneously reproducing the prope rties of local 
spirals (jPi ontck fc Stcinmctz 2011:_Governa to et al.l [20091 : 
IScannapieco et al.ll2009l : iBrooks et al]i2009bl ') . Moreover, re- 
cent observational studies indicate that star formation rate 
densities in local galaxies as well as at high redshift corre- 
late linearly with the surface density of the molec ular gas, 
with no evolution in the molecular SK relatio n (Bi giel et al.l 
I2OO8I : iDaddi et al.] I2OI0I : iGenzel et ai]|2010f) . This is also 
true for galaxies with very clumpy star formation, as ex- 
pected for a linear dependence with gas density. Only in- 
teracting galaxies undergoing a signiflcant starburst seem 
to show an increased sta r formation efficiency (Dad di et al.. 
I2OIOI : iGenzel et al.|[201ol 'l. Neither the SAMs nor the simu- 
lations presented here include these effects, although we are 
working on updated models that will do so. 

A second major difference between the SAMs and the 
simulations is in the mode in which galaxies acquire most of 
their stellar mass. We distinguish between "in situ" growth, 
due to stars that form out of cold gas within the galaxy 
in question, and "accretion" of stars that formed in exter- 
nal galaxies and are accreted via mergers. This distinction 
is important because it may dete rmine the ch aracteristic 
size evolution of early type galaxies ( Khochfar fc Silk 20061: 



20081: lOser et all I2OI0I : iFeldmann eT^ I2OI0I : IZehavi et all 



2011f ). Although the ratio of in situ to accreted stars always 



iNaab et al.|[2009bl : IGuo et al.|[2oTll : 



Covington et al.ll201 



It has been shown previously that massive galaxies in the 
simulations and semi-analytic models have an early phase of 
growth dominated by the in situ mode, and then switch over 
to a mode that is dominated by growth through accretion 
jPe Lucia et al.ll2006l : iKhochfar fc sii^l2006cl : IGuo fc Whit3 



decreases with time in the SAMs, in qualitative agreement 
with the simulations, in the SAMs in situ growth dominates 
over accretion in halos of all masses at all times. Examining 
the absolute mass in stars formed in situ or accreted, we find 
that this is primarily because the SAMs predict much less 
mass growth through accretion — the in situ mass evolu- 
tion agrees fairly well with the simulation results. The com- 
parison between SAMs with different physical processes in- 
cluded gives us further insights into the origin of the discrep- 
ancy: 1) Increasing the efliciency of SN FB (SN and SNWM 
SAMs) reduces the accreted mass because star formation 
is suppressed in the low-mass satellites that eventually get 
accreted. 2) Including radio mode AGN FB (FULL SAM) 
reduces the in situ growth in massive galaxies at late times, 
because it shuts off the fuel supply for in situ star formation 
in these objects. Interestingly, increasing the SF efliciency 
in the SAM does not affect the fraction of in situ versus 
accreted mass, presumably because all galaxies are boosted 
equally. 

In summary, it is likely that the contributions from in 
situ and accreted stars are currently incorrectly predicted in 
both the simulations and SAMs, for the following reasons. In 
order to match observations, the simulations presented here 
clearly require both a process that suppresses star formation 
in low mass objects at all redshifts (such as SN-driven winds) 
and one that can shut off residual cooling and quench star 
formation in massive galaxies at late times (such as radio 
mode AGN FB). The former will reduce the accreted mass, 
while the latter will decrease the in situ mass in massive ob- 
jects at late times. In addition, if the SFE is higher in high 
redshift galaxies than at late times, this is likely to increase 
the accreted mass. On the other hand, the SAMs presented 
here, like many SAMs in the literature, make the assump- 
tion that hot gas from the halo can only be accreted onto 
the central galaxy. This rapidly truncates the star forma- 
tion in satellite galaxies and is likely to artificially decrease 
the accreted mass fraction. As a possible solution, strangu- 
lation of satellite galaxies might be delayed if the depth of 
the potential well of the sub-halo is deep enough to retain 
the hot halo gas for a longer time (e.g. Khochfar fc Ostriker 
2008). In this case, the accreted objects may have longer on- 
going cooling resulting in longer on-going star formation. In 
addition, the SAMs neglect gravitational heating, which is 
clearly important in the simulations and should reduce the in 
si tu growth in ma s sive h alos at late times. The recent study 
of iFontanot et all (120111 ) has shown that the S08 SAM and 
other similar SAMs in which radio mode feedback is used 
to solve the overcooling problem overpredict the fraction 
of radio loud galaxies compared with observations, and the 
implemented dependence of radio luminosity on stellar and 
halo mass is too steep. Gravitational heating could play a 
similar role jKhochfar fc Ostrikei|[20oi : iDekel fc BirnboimI 
,20081 : iBirnboim fc Dekej|2010l ') and thereby reduce the need 
for such strong radio mode heating. 

A third important result is that the cooling recipe im- 
plemented in these SAMs (which is widely used in many 
SAMs) overpredicts the overall accretion rate of gas by a 
factor of 1.5 (in low mass halos) to four (in massive ha- 
los). If the accretion rates in the simulations are accurate, 
this implies that the SAMs that match present day galaxy 
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properties are compensating for this "extra" accretion by 
artificially "tuning up" the feedback. Moreover, when we di- 
vide the gas accretion int o "hot mode" and "cold mode" 
iBirnboim fc Dekel 2003a : Keres et aLlbOOSbl ). we find that 
the SAMs systematically overestimate the hot mode growth 
(by up to an order of magnitude) and underestimate the cold 
mode. As well, in the SAMs without metal cooling, the cold 
mode shuts off completely at low redshifts (the shutoff red- 
shift depends on halo mass), while in the simulations it de- 
clines smoothly but continues to low redshifts. This is likely 
to be a result of the fact that, in the SAMs, the criterion for 
discriminating between hot and cold mode is based on the 
assumption of smooth, spherical halos, while in simulations 
cold gas can stream into halos along cold, dense filaments. 
In addition, in the SAMs, gas is assumed to accrete either in 
cold mode or hot mode, but simultaneous cold and hot mode 
accretion is not allowed. In the simulations, the dense, cold 
streams can penetrate deep into the diffuse hot halos, allow- 
ing for both accretion modes to occur within the same halo 
l|Keres et al.ll2005bl , l2009l : lBrooks et al.ll2009al ). We consider 
this discrepancy to be a general feature of currently exist- 
ing SAMs (see section and for example, similar results 
were found in the study of Lu et al. (2010). In order to over- 
come the weakness of SAMs in this respect, they proposed 
a new model that explicitly incorporates cold-mode accre- 
tion independent of the hot halo gas. By fitting the hot and 
cold gas fraction in simulations as a function of redshift and 
halo mass and assuming accretion onto the galaxy within a 
free-fall time, they calculate the accretion rate of the cold 
component and thus, achieve a better match of their SAMs 
to simulations. However, motivated by hydrodynamical sim- 
ulations, implementing gravitational heating might provide 
a more promising mechanism as it treats gas inflow and heat- 
ing in a self-consistent way. Following Khochfar & Ostriker 
(2008), the heating of the intracluster medium (ICM) would 
be calculated by the net surplus of gravitational potential 
energy released from gas that has been stripped from in- 
falling satellites. Gravitational heating is found to be an ef- 
ficient heating source for massive dark matter halos, where 
it prevents cooling, and becomes especially important at late 
times. Therefore, if we generally assumed cold infalling gas 
and heating of gas due to the binding energy of the halo 
potential, this might imply an automatic, less efficient heat- 
ing at high z and for low mass halos. This should result 
naturally in larger cold accretion fractions for these objects 
and thus, also in a better match of the hot gas content for 
high-redshift galaxies. 

We compared both the simulations and SAMs to two 
key observational constraints at z = 0, 1, and 2: the re- 
lationship between dark matter halo mass and stellar mass, 
and the relationship between stellar mass and star formation 
rate. We found that the simulations predict stellar masses 
that are too large for their halo masses at all redshifts. The 
stellar masses are too high by a factor of a few for massive 
halos ( ^ 10^^ Mq), and by an order of magnitude or more 
for lower mass halos. The SAM results for the NF model are 
qualitatively similar, although the stellar masses at high red- 
shift are lower, due to the lower star formation efficiencies, 
as already discussed. In the SAMs, the curvature in the em- 
pirical Mgai — Afhaio relation can be achieved by including 
supernova driven winds, which suppress star formation in 
low mass halos, and radio mode AGN FB, which suppresses 



star formation in high mass halos. We found that the SFR 
at a given stellar mass were too low in the simulations at 
all redshifts z ^ 2, probably because of the overly efficient 
SF at higher redshifts and the resulting gas consumption. 
SF galaxies in all SAMs were found to lie on the same re- 
lation, with supernova feedback shifting the galaxies along 
the relation. The SAMs (both FULL and NF) showed better 
agreement with the observed rhstar — Afgai relation than the 
simulations, but have SFR at a given stellar mass that are 
about a factor of ~ 2 too low at high redshifts. 

In final summary, we conclude that on the one hand, 
we are encouraged by the robustness of SAMs as a tool 
for exploring the qualitative effects of varying the phys- 
ical ingredients of galaxy formation. On the other hand, 
we have also identified several important areas where the 
quantitative accuracy of fundamental physical recipes in 
the SAMs should be improved, and several physical pro- 
cesses that are missing in our SAM but which should be 
included. Additionally, there is a tendency to treat numer- 
ical simulations as "truth", but we have shown that key 
predictions of these simulations are sensitive to uncertain 
sub-grid recipes. We have also highlighted several physi- 
cal processes that are neglected in the simulations studied 
here, but which appear to be crucial in order to under- 
stand the properties of real galaxies. These include more ef- 
fective implementation of supernova-driven winds, chemical 
enrichment and metal cooling, and a self-consistent treat- 
ment of the growth of and feedback from black holes. Of 
course, these are hardly new suggestions, and considerable 
progress has been made recently in all of th e se areas (e.g. 
Pi Matteo et al.l l2005l: IC attaneo et al.l liooi : 'Siiac ki et al] 
20071: lOppenheimer fc Dave .20081: iBooth fc Schavd l2009l: 



Scannapieco et al .1120091: 1 iGovernato et al .Il2010l:l Schave et al.l 



2010bl : ISawala et allboiol : lOstriker et al]|2010l '). In conclu- 
sion, we suggest that using these two complementary tech- 
niques (SAMs and hydrodynamic simulations) together in 
close coordination may provide the most powerful approach 
to understanding galaxy formation and evolution for the 
near future. 
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Figure Al. Upper panel: comparison of aggregation history for 
halo M0501, with total mass 3.0 X IO^^A/q, for two different reso- 
lutions (2x: red dotted line and 4x: black solid line). Intermediate 
panel: Number of minor mergers (> 10 : 1) as a function of look- 
back time for the two different resolutions. Lower panel: Number 
of major mergers (< 10 : 1) as a function of lookback time for the 
two different resolutions. 



APPENDIX A: EFFECTS OF NUMERICAL 
RESOLUTION 

Al Dark matter component 

For one high and low mass halo with final masses of Mhaio = 
3 X 10^^ Mq (M0501) and of A4aio = 8 x lO^M© (M1646), 
we have performed re-simulations with 4x the original spa- 
tial resolution. We traced back the particles that are closer 
than 2 x r2oo to the center of the halo in any of our snap- 
shots and replace them with dark matter and gas particles 
of higher resolution, achieving a 16x better mass resolu- 
tion in the high resolution region than in the original sim- 
ulation: moM ~ 2.5 X 10^Mq/i~^ and mcas ~ m-star ~ 
5 X 10^Mq/i-\ The mass aggregation history of the dark 
matter component in both halos is very similar for the two 
different resolution limits as shown in the upper panel of 
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Figure A2. Same as Fig. lAll but for halo M1646, with a smaller 
mass of 8.0 X IO^Mq. 



Figs. IAll and lA2l In the middle and lower panel of Figs. lAll 
and lA2l one can see from the histograms of major (> 1 : 10) 
and minor mergers (< 1 : 10) that the number of major 
mergers stays almost the same, while as expected, we can 
resolve more minor mergers (< 1 : 10) in the higher res- 
olution case. This suggests that, since galajcy formation in 
the SAMs is mainly influenced by the accretion history of 
the main halo and by major merger events (> 1 : 10), the 
results from the SAMs are expected to be well-converged in 
our re-simulations. 
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Figure A3. Evolution of the baryonic components for two dif- 
ferent mass resolutions (2x: dotted lines, 4x: solid lines) for halo 
M0501. The upper panel shows the evolution of the stellar compo- 
nent, the middle panel the evolution of the cold gas fraction and 
the lower panel shows the hot halo gas component. Black lines 
illustrate the re-simulations, and colored lines show the different 
SAM versions (see main text). 



in both the SAMs and the simulations. Therefore, due to 
computational costs, we restrict our main study to 2x re- 
simulations. 



A2 Baryonic components 

Figs. IA3I and IA4I show explicitly the evolution of the star, 
cold and hot gas mass in SAMs and simulations based on 
the trees from the two re-simulated individual halos for 2 x 
and 4 x higher resolution. As expected from the similar 
evolution of the dark matter component, for both the low 
and the high mass halo re-simulations, we find no significant 
difference in the evolution of the central galaxy/main halo 
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Figure A4. Same as Fig.|A3l but for halo M1646. 



